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Abstract 



We investigate the associated production of heavy quarkonia, with angular- 
momentum quantum numbers 2S+1 Lj = %, 3 S"i, x Pi*Pj (J = 0, 1, 2), and photons, 
Z bosons, and W bosons in photon-photon, photon-hadron, and hadron-hadron 
collisions within the factorization formalism of nonrelativistic quantum chromody- 
namics providing all contributing partonic cross sections in analytic form. In the 
case of photoproduction, we also include the resolved-photon contributions. We 
present numerical results for the processes involving J/tj} and \cj mesons appropri- 
ate for the Fermilab Tevatron, CERN LHC, DESY TESLA, operated in the e+e" 
and 77 modes, and DESY THERA. 
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1 Introduction 



Since its discovery in 1974, the J /if) meson has provided a useful laboratory for quantita- 
tive tests of quantum chromodynamics (QCD) and, in particular, of the interplay of per- 
turbative and nonperturbative phenomena. The factorization formalism of nonrelativistic 
QCD (NRQCD) fl|] provides a rigorous theoretical framework for the description of heavy- 
quarkonium production and decay. This formalism implies a separation of short- distance 
coefficients, which can be calculated perturbatively as expansions in the strong-coupling 
constant a s , from long-distance matrix elements (MEs), which must be extracted from 
experiment. The relative importance of the latter can be estimated by means of velocity 
scaling rules, i.e., the MEs are predicted to scale with a definite power of the heavy-quark 
(Q) velocity v in the limit v <C 1. In this way, the theoretical predictions are organized 
as double expansions in a s and v. A crucial feature of this formalism is that it takes into 
account the complete structure of the QQ Fock space, which is spanned by the states 
n = 2S + l l}j> with definite spin S, orbital angular momentum L, total angular momentum 
J, and color multiplicity c = 1,8. In particular, this formalism predicts the existence of 
color-octet (CO) processes in nature. This means that QQ pairs are produced at short 
distances in CO states and subsequently evolve into physical, color-singlet (CS) quarkonia 
by the nonperturbative emission of soft gluons. In the limit v — > 0, the traditional CS 
model (CSM) is recovered. The greatest triumph of this formalism was that it was 
able to correctly describe the cross section of inclusive charmonium hadroproduction 
measured in pp collisions at the Fermilab Tevatron || , which had turned out to be more 
than one order of magnitude in excess of the theoretical prediction based on the CSM. 

In order to convincingly establish the phenomenological significance of the CO pro- 
cesses, it is indispensable to identify them in other kinds of high-energy experiments as 
well. Studies of charmonium production in ep photoproduction, ep and uN deep-inelastic 
scattering, e + e~ annihilation, 77 collisions, and 6-hadron decays may be found in the 
literature; see Ref. || and references cited therein. Furthermore, the polarization of char- 
monium, which also provides a sensitive probe of CO processes, was investigated ||[7J . 
Until very recently, none of these studies was able to prove or disprove the NRQCD fac- 
torization hypothesis PJ. However, preliminary data of 77 — > J /if) + X taken by the 
DELPHI Collaboration || at LEP2 provide first independent evidence for it |J. 

In most NRQCD studies of charmonium production in high-energy particle collisions, 
the charmonium state is produced either singly or in association with a hadron jet (j), 
which originates from an outgoing quark (g) or gluon (<?), because this is bound to yield 
the largest cross sections 0,[|[|[7|||. However, also the associated production of J /if) 
mesons and prompt photons (7) was considered, in photon-photon f]10| ,|TT1 , pL2|, pL3[| , photon- 
hadron fl^ , [15|j , and hadron-hadron collisions p~5|JT6|] . Although the cross sections of the 
latter processes are suppressed by a factor of a/a s , where a is Sommerfeld's fine-structure 
constant, relative to the case of J/ip+j associated production, these processes exhibit some 
attractive features. On the one hand, J /if) + 7 associated production in photon-photon 
collisions dominantly proceeds through the CS partonic subprocesses of direct and doubly- 
resolved photoproduction, 77 — > cc 3 Si 7 and gg — > cc 7, respectively, and thus 



2 



allows for an independent determination of the CS ME ((9 J// ^[ 3 Sj 1 ' ) ]), which is usually 
extracted from the measured leptonic annihilation rate of the J/ip meson. On the other 
hand, J/tp + 7 associated production in photon-hadron collisions provides a good handle 
on the CO mechanism because the partonic subprocesses of direct photoproduction, jg — > 
cc[n]7 with n = ^ , 3 S{ , 3 P j , are pure CO processes [[T3],[nj. Finally, for small J/ip 
transverse momentum, J/ip + 7 associated hadropro duct ion chiefly proceeds through the 
CS partonic subprocess gg — > cc 3 Si 7 and thus lends itself as a clean probe of the 
gluon density inside the proton [|l5|,|n|. 

The purpose of this paper is to study the associated production of heavy quarkonia 
with the intermediate bosons, W and Z, in photon-photon, photon-hadron, and hadron- 
hadron collisions to leading order (LO) in the NRQCD factorization formalism. We 
consider all experimentally established heavy quarkonia, which are classified by their 
angular-momentum quantum numbers 2S+1 Lj = l So, 3 Si, 1 Pi, 3 Pj with J = 0,1,2. In 
the case of charmonium, these correspond to the r] c , J/ip (ip 1 , ip(3S), etc.), h c , and Xcj 
mesons, respectively. We present all contributing partonic cross sections in analytic form. 
As a by-product of our analysis, we recover the corresponding formulas for the associ- 
ated production of heavy quarkonia with prompt photons, by appropriately adjusting 
the couplings and mass of the Z boson. In our numerical study, we concentrate on those 
charmonia which can be most straightforwardly identified experimentally, namely the J / ip 
and Xcj mesons, through their leptonic and radiative decays, respectively. Specifically, we 
assess the feasibilities of the pp synchrotron Tevatron, which is in operation at Fermilab, 
the pp synchrotron LHC, which is under construction at CERN, the e + e~ linear collider 
(LC) TESLA, which is being developed and planned at DESY and can be operated in 



the e + e [17]] and 77 ]18| modes, and the pe collider THERA, which uses the TESLA 



lepton beam and the HERA proton beam [I9|, to produce J/ip and Xcj mesons in asso- 



ciation with prompt photons, W bosons, and Z bosons. These colliders all operate with 
high luminosities at high energies and, at first sight, have the potential to produce these 
final states with reasonable rates. Furthermore, if one selects the leptonic decays of the 
J/ip meson and intermediate bosons, then these final states lead to spectacular signals, 
consisting of energetic isolated prompt photons, charged leptons, and possibly missing 
four-momentum, which should be easy to detect even at hadron colliders. A first step in 



this direction was recently undertaken in Ref. ||20|| , which contains a numerical study of 
T + W and T + Z hadroproduction at the Tevatron and the LHC. 

For completeness, we also provide numerical results for the cross sections of the pro- 
duction of J/ip and Xcj mesons together with standard-model (SM) Higgs bosons, which 
are suppressed by the smallness of the charm Yukawa coupling. In fact, explicit analysis 
reveals that these cross sections are too small to produce any signals for the envisaged 
luminosities of the considered colliders, even in the case of the LHC. By the same token, 
the observation of events where a scalar boson is produced together with charmonium 
would signal new physics beyond the SM. 

This paper is organized as follows. In Sec. |2], we present our analytic results and 
explain how to evaluate the cross sections of the various processes enumerated above. 
The contributing partonic cross sections are collected in the Appendix. In Sec. ^ we 
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present our numerical results. Our conclusions are summarized in Sec. [|. 



2 Analytic results 

In this section, we explain how to calculate the cross sections of the associated production 
of a heavy quarkonium state C and an electroweak boson D by photoproduction in e + e~ 
and pe ± collisions and by hadroproduction in pp and pp collisions. 

To start with, we consider e + e~ collisions, where photons are unavoidably generated by 
hard initial-state bremsstrahlung. At a high-energy e + e~ LC, an additional source of hard 
photons is provided by beamstrahlung, the synchrotron radiation emitted by one of the 
colliding bunches in the field of the opposite bunch. Bremsstrahlung and beamstrahlung 
both occur in the e + e~ mode, and their energy spectra and must be added coherently. 
The highest possible photon energies with large enough luminosity may be achieved by 
converting the e + e~ LC into a 77 collider via back-scattering of high-energetic laser light 
off the electron and positron beams. 

The photons can interact either directly with the quarks participating in the hard- 
scattering process (direct photoproduction) or via their quark and gluon content (re- 
solved photoproduction). Thus, the inclusive process e + e~ — > e + e~CD + X, where X 
denotes the hadronic remnant, receives contributions from the direct, singly-resolved, and 
doubly-resolved channels. All three contributions are formally of the same order in the 
perturbative expansion and must be included. This may be understood by observing that 
the parton density functions (PDFs) of the photon have a leading behavior proportional 
to a ln(M 2 /AQ CD ) oc a/a s , where M is the factorization scale and Aq CD is the asymptotic 
scale parameter of QCD. 

Let us now describe how to calculate the cross section. We take the electron and 
positron to be massless and denote the masses of C and D by mc and mo, respectively. 
Let be the energy of the initial state, yc and yn the rapidities of C and D, respectively, 
and Pt their common transverse momentum in the center-of-mass (cm.) frame of the 
collision. Invoking the Weizsacker- Williams approximation (WWA) for electromagnetic 



bremsstrahlung [21], its analogues for beamstrahlung [E3,Ej] and Compton back-scattering 



25,26], and the factorization theorems of the QCD parton model p2[ and NRQCD p 



the differential cross section of e + e — > e + e CD + X can be written as 

d 3 a , , 



dp\ dy c dy D 



e e 



e + e CD + X) = / dx + f 1 / e {x + ) \ dx- / 7 / e (x_) 



X 



da 



E Xafa/ y (x a , M)x b f b/l {x b , M) (O [n]) — (ab -> QQ[n] + D) 



a.b.n 



where it is summed over the active partons a,b — J,g,q,q, / 7 / e (^±) are the photon flux 
functions, / a / 7 (x a , M) and f b/l (x b , M) are the PDFs of the photon, {O c [n}) are the MEs 
of C, (da/dt)(ab — > QQ[n] + D) are the differential partonic cross sections, and x± and 
x a ,b are the fractions of longitudinal momentum that the emitted particles receive from 
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the emitting ones. With the definition / 7 / 7 (x 7 , M) = 8(1 — x 7 ), the direct, singly-resolved, 
and doubly-resolved channels are all accommodated. We have 

m% exp(±y c ) + m% exp(±y D ) 



x± 



%a,b 



s 

in? exp(±?/c*) + exp(±yc 



x± 



(3) 



where = \jvn? c + p\ is the transverse mass of C and similarly for D. 

The partonic Mandelstam variables s = (p a +Pb) 2 , t = (p a — pc) 2 , and u = (p a 
can be expressed in terms of pr, yc, and yo, as 



Pd) 



s 

t 
u 



(m 



t) 



+ (m^J + 2m^my coshf^c — y D ), 



-p 2 T - m%m^ exp(y D - y c ), 
-p 2 T - m%m^ exp(y c -yo), 



respectively. Notice that s + t + u = m 2 c + m 2 D and sp\ 
allowed ranges of S, px, yc, and yr> are 



tu ■ 



i i 
m c m D . 



(4) 

The kinematically 




<p T < 



\yc\ < Arcosh 



S + m 2 r 



c 



m 



D 



2VSm c 



T 



hi 



Vs- 



171% exp( 



-yc) 



rn 



D 



< Vd < hi 



\f~S~ — exp(y, 



c 



where A(x, y, z) = x 2 + y 2 + z 2 — 2(xy + yz + zx) is Kallen's function. 

Sometimes it may be interesting to consider the cross section differential in the CD 
invariant mass tticd — \fs rather than in p?- Using 

r cosh(?/c - y D ) - m 2 CD 



D 



(5) 



Pt 



2sinh 2 (y c - y D ) 



(6) 



where r 



m C D + ( m c 



2 D ) 2 smh 2 (y 



c 



y D ), we find 



dp 2 T 
dm 2 CD 



Equations (Bl) and 



m 2 CD cosh(y c -y D )-r 
2rsinh 2 (y c - Vd) 
are regular for y c = y D and read then 
A (S, m 2 c ,m 2 D ) 



(7) 



Pt 



AS 



dp 



dm 2 CD 



m 



c 



m 



rn 



CD 



5 



respectively. Multiplying Eq. (|1|) with Eq. (|7|), we obtain (d 3 a jdrv? CD dye dyn) (e + e~ — > 
e + e~CD + X). The kinematically allowed range of mcD is trie + rri£> < tticd < y/S- 

We work in the fixed-flavor-number scheme, i.e., we have rif = 3 active quark flavors 
q = u, d, s if Q — c and rif = 4 active quark flavors q = u,d, s,c if Q = b. As required 
by parton-model kinematics, we treat the q quarks as massless. To LO in v, we need to 
include the cc Fock states n = ^ , tfj , 3 S{ 8 \ l P? ] if H = Vc ; n = 3 S[ 1) , 1 S^\ 3 S[ 8 \ 3 Pj 8) if 
H = J/ij>, 4>(3S), . . .; n = ^ 8) if H = h c ; and n = 3 Pj 1} , 3 S[ 8) if H = x cJ , where 
J = 0, 1, 2 [ED]. Their MEs satisfy the multiplicity relations 



QipinS) 
QXcJ 
QXcj 



3p(8) 
3p(l) 
'35.(8) 



(2J + 1) (o^ nS) 



3p(8) 



(2 J + 1) (O Xc0 
(2J + 1) (C xc0 



3p(l) 
M) 

3 5{ 8) ' 



(9) 



which follow to LO in v from heavy-quark spin symmetry. The assignments for the various 
bottomonia are analogous. 

The cross section of pe^ 1 — ► e ± CD + X in photoproduction emerges from Eq. ([I]) 
by substituting f 1 / e {x + ) = 5(1 — x + ) and replacing the photon PDFs f a /-y{x a , M) with 
their proton counterparts f a / p (x a , M). It receives contributions from the direct and re- 
solved channels. Here, it is understood that / 7 / e (x_) only accounts for electromagnetic 
bremsstrahlung. To conform with HERA conventions, we take the rapidities to be posi- 
tive in the proton flight direction. The rapidity y l ^ h of C in the laboratory frame, where 
the proton and electron have energies E p and E e , respectively, is related to yc by 



y l S h 



!, E p 
yc+ 2 ln E-/ 



(10) 



E p E e 



and similarly for D. The cm. energy is \/~S = 2 1 

In a second step, the cross section of pp — > CD+X in hadroproduction is obtained from 
the one of pe ± — > e ± CD + X in photoproduction by substituting / 7 / e (x_) = 8(1 — xJ) and 
replacing the photon PDFs /6/ 7 (x&, M) with their antiproton counterparts fb/p(xb, M). In 
hadron-collider experiments, the pp cm. frame and the laboratory one usually coincide. 
The cross section of pp — > CD + X in hadroproduction is accordingly evaluated with 
f b/p (x b ,M). 

We now turn to the partonic subprocesses ab — > QQ[n]D. The differential cross 
section of such a process is calculated from the pertaining transition-matrix element T 
as da/dt = \T\ 2 / (167rs 2 ), where the average is over the spin and color degrees of freedom 
of a and b and the spin of D is summed over. We apply the covariant-projector method 
of Ref. P^J to implement the QQ Fock states n according to the NRQCD factorization 
formalism [[J. 

The following partonic subprocesses contribute to LO in a s and v: 



77 - QQk {1) }N, 

ig^QQ[^ 8) }N, 



(11) 
(12) 
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gg^QQ[qU]N, (13) 

gg^QQ[q^]N, (14) 

qq^QQ[^]N, (15) 

quq d ^QQ[ 3 S{ S) ]W + , (16) 

-> gg fe{ 8) W~ (17) 



where q — u,d, s, c; g u = w, c; qd = d, s; N = 7, Z; and <j = ^o, 3 £i, ^1, 3 Pj with J = 0, 1, 2. 
For the reason explained above, q = c and q u = c must not be included if Q = c. The 
processes 77 — > QQ[q^]N and jg — ► QQ[?W]iV are forbidden by color conservation. 
Furthermore, the processes gg — > QQfc are prohibited because the g-quark line is 
connected with the g-quark line by one gluon, which transmits color to the QQ pair. For a 
similar reason and due to the fact that the W boson must be emitted from the initial-state 
quarks if a QQ pair is to be produced, processes ( |16|) and fli~7|) only come with n = 3 S[ 8 \ 

In the Higgs-boson case D = H, the contributing partonic subprocesses are analogous 
to Eqs. (|llf)-(|l5"[), except that, due to charge- conjugation invariance, 77 — > QQ[q^]H, 
19 -> QQk {8) )H, and gg -> QQ^^iZ" are forbidden for ? = 3 Si, iPi, and gg -> QQ[^]H 
is forbidden for ^ = ^o, 3 Pj. 

The differential cross sections da/dt of processes fllTD-fllTD are listed in the Appendix 
We combine the results proportional to the CO MEs /C^(™ 5 ) 
exploiting the multiplicity relations of Eq. (|9[). We include large logarithmic corrections 
due to the running of the fine-structure constant from the Thomson limit to the elec- 
troweak scale by expressing the gauge couplings in terms of Fermi's constant Gp, as 
g = 2 1 / 4 G 1 J 2 m z and g' = 2 3 / 4 G 1 /' 'mw- Furthermore, we define the Zqq vector and 
axial- vector couplings as v q = I 3 — 2e q sm 2 9 w and a q = I 3 , respectively, where e q is the 
fractional electric charge of quark q, I 3 is the third component of weak isospin of its left- 
handed component, and 6 W is the weak mixing angle, which is fixed by cos 9 W = m w /m z . 
The results for processes (Tf2] ) and flT3"|) may be obtained from those for process (JTTJ) by 
adjusting the overall color factors as specified in the Appendix. The results for N = 7 
emerge from the ones for iV = Z by adjusting the couplings and mass of the Z boson, by 
substituting g = e, v q = e q , a q = 0, and m z = 0, where e = y/Aita. is the electron charge 
magnitude, and they are not presented separately. They vanish for processes (pTT|) — (|1~3|) 
with q = l So, 3 Pj and for process (p~5|) with q = X P\ because the photon has no axial- vector 
coupling. 

In the case of N = 7, the differential cross sections da/dt of processes (|TT|) — ( with 



q = 1 P 1 and of process (|l^) with q — 3 5i, 3 Pj are plagued by infrared singularities in the 
limit of the final-state photon being soft, i.e., for t — u — 0. In addition, process (|T5|) 
with q = 3 S\ also suffers from a u- or t-channel singularity if the final-state photon is 
hard and collinear to the initial-state q or q quarks, i.e., for t < u = or u < t = 0, 
respectively. Owing to the identity ])\ = tu/s, both the soft and collinear limits entail 
p T — > 0. By the same token, in the evaluations of Eq. ([!]) and its counterparts for lepton- 
hadron and hadron-hadron scattering, all these singularities can be avoided by imposing 
a lower cut-off on p?. On the other hand, detailed inspection of these equations reveals 
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that da/dpT is finite in the limit px — > if da/dt is.[] The cases D = Z,W, H are devoid 
of such singularities because m D acts srulator. 

To our knowledge, the formulas presented in the Appendix cannot be found else- 
where in the literature. However, the literature contains analytic results for the partonic 
subprocesses pertinent to the case of C = J/ip and D — 7, namely, for processes (|TT| ) 
and (p~3|) ||T4| , |16| with q = 3 Si, and for processes (|T4|) and 
[114, 1611 . We agree with them, except for Eq. (7) of Ref. JTT] 



H3, 

1 C 3 c 
*->0j >->l) 



(up m 



related comments in Refs. P, P^JT3[ ]) and the term proportional to (0% ( 3 Po)/ i n Eq. (7) 
of Ref. H-Q 



T§ with 
(see also 



3 Numerical results 



We are now in a position to present our numerical results. We focus our attention on 
the cases C = J/ip, Xcj- These charmonia can be efficiently identified experimentally, and 
their MEs are relatively well constrained J7|. The predicted cross-section distributions 
for the ip' mesons are similar to those for the J/ip mesons, but their normalization is 
somewhat suppressed due to smaller MEs 0. The r) c meson is more difficult to detect 



experimentally, and the h c meson is poorly known |28 



We first describe our theoretical input and the kinematic conditions. We use m c = 
m c /2 = 1.5 GeV, m w = 80.423 GeV, m z = 91.1876 GeV, G F = 1.16639 x 1(T 5 GeV~ 2 , 
a = 1/137.036, and the LO formula for ai n/ ^(/i) with rif = 3 active quark flavors f28fl . 
We assume ttih = 115 GeV, a value just above the 95%-confidence-level lower bound 
on win from direct Higgs-boson searches and in agreement with tantalizing hints for 
the direct observation of the Higgs-boson signal delivered towards the end of the LEP2 
running phase. As for the photon PDFs, we use the LO set from Gliick, Reya, and 
Schienbein (GRS) p9|| , which is the only available one that is implemented in the fixed- 
flavor-number scheme, with rif = 3. As for the proton PDFs, we use the LO set from 

Martin, Roberts, Stirling, and Thorne (MRST98LO) Q. We adopt A§£ D = 204 MeV 

from Ref. which happens to precisely correspond to Aq^ d = 174 MeV, the value 
employed in Ref. [30||, if the matching scale is taken to be m c . We choose the renormal- 

ization and factorization scales to be [i = M = if D = 7 and to be ft = M = ^Jm^m^ 
if D = Z,W,H. As for the J/ip and Xcj MEs, we adopt the set determined in Ref. 0] 
using the MRST98LO proton PDFs. Specifically, (p J ^ [ 3 S{ 1] ]) and (O* c0 [ 3 P (1) ] ) were 
extracted from the measured partial decay widths of J/ij) — > l + l~ and Xc2 — > 77 p8| , 



respectively, while lO J ^ l S$ 



7(8) 



3o(8) 



qJ/4> 



3 Pr 



(8) - 




and (O Xc0 



35.(8) 



were 



fitted to the transverse-momentum distributions of J / ip and XcJ inclusive hadroproduction 

1 In this sense, we disagree with the statement made in Ref. that process ([u]) with = ^o, 3 -Po, 3 P2 
can produce collinear or infrared singularities. The triple-gluon vertex present in the contributing Feyn- 
man diagrams is innocuous because the virtual gluon has a timelike virtuality in excess of vr? c . 

2 There are two obvious typographical errors in Eq. (8) of Ref. [Q: (2m;?) and the second appearance 

Of (O; 



Si) ) should be replaced by (2m c ) and 



qJJiI> respectively. 
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and the cross-section ratio cr Xc2 /a Xcl |31] measured at the Tevatron. The fit results for 



Qj/lp 



lo(8) 
a 



and (O j ^ 



3 P t 



(8) 



are strongly correlated, so that the linear combination 



\ i j / m j \ i j / 



with a suitable value of r 
linear combination of 



is quoted. Unfortunately, Eq. (|TD is sensitive to different 



and (O j ^ 



3 Pr 



than appears in Eq. (|i8|) . In 



want of more specific information, we thus make the democratic choice (O j ^ 
(r/m 2 ) (OW 



3 P< 



(8) 



lc(8) 
D 



= MW/2. 

We now discuss the photon flux functions that enter our predictions for photopro- 
duction at TESLA and THERA. The energy spectrum of the bremsstrahlung photons is 
well described in the WWA |ZT| by Eq. (27) of Ref. |J5). We assume that the scattered 
electrons and positrons will be antitagged, as was usually the case at LEP2, and take the 
maximum scattering angle to be # max = 25 mrad ||33|| . The energy spectrum of the beam- 
strahlung photons is approximately described by Eq. (2.14) of Ref. ||24|| . It is controlled 
by the effective beamstrahlung parameter T, which is given by Eq. (2.10) of that refer- 
ence. Inserting the relevant TESLA parameters for the \fS = 500 GeV baseline design 



specified in Table 1.3.1 of Ref. |l7j in that formula, we obtain T = 0.053. In the case 
of the e + e~ mode of TESLA, we coherently superimpose the WWA and beamstrahlung 
spectra, while, in the case of THERA, we only use the WWA spectrum. Finally, in the 
case of the 77 mode of TESLA, the energy spectrum of the back-scattered laser photons 

It depends on the parameter k 



is given by Eq. (6a) of Ref. 25 . It depends on the parameter k = s ei /m A e — 1, where 
y/s^ is the cm. energy of the charged lepton and the laser photon, and it extends up to 
i max = k/(k + 1), where x is the energy of the back-scattered photons in units of y/~S/2. 
The optimal value of k is k = 2 (l + \/2j ~ 4.83 ||26|| , which we adopt; for larger values 
of k, e + e~ pairs would be created in the collisions of laser and back-scattered photons. 

Since our study is at an exploratory level, we refrain from presenting a quantitative 
estimate of the theoretical uncertainties in our predictions. However, experience from 
previous analyses of charmonium production within the NRQCD factorization formalism 



34}| leads us to expect relative errors of the order of ±50%. 



Table [I] gives the cm. energies Vo, design luminosities L, numbers of dedicated ex- 
periments, and total cross sections o\ corresponding to a yield of one signal event per 
year of operation for the various colliders considered. At THERA, the proton and lepton 
energies in the laboratory frame are planned to be E p = 1 TeV and E e = 250 GeV fl9|| . 
One year of operation is usually taken to be 10 7 s, which corresponds to a duty factor of 
approximately 33% and can be achieved, e.g., by 200 days of running with an efficiency 
of 60%. Thus, the figures for L in units of 10 32 cm~ 2 s -1 presented in Table [I] equally 
correspond to the integrated luminosity J dt L in units of fb _1 per year and experiment. 
At the Tevatron and the LHC, the measurements can be performed simultaneously by 
two dedicated experiments. We thus obtain the values of o\ specified in Table [TJ 

We are now in a position to present our numerical results. Figures |l]-|5| are devoted 
to e + e~ — ► e + e~CD + X in the e + e~ and 77 modes of TESLA, to pe^ 1 — - > e ± CD + X at 
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Table 1: Cm. energies v5, design luminosities L, numbers of dedicated experiments, and 
total cross sections U\ corresponding to a yield of one signal event per year of operation 
for the considered colliders. 



Collider 


[TeV] 


L [10 32 cm-V 1 ] 


No. of experiments 


<Ti [fb] 


TESLA e + e~ mode 


0.5 


340 


1 


0.0029 


TESLA 77 mode 


0.5 


60 


1 


0.017 


THERA 


1 


0.041 


1 


24 


Tevatron Run II 


2 


2 


2 


0.25 


LHC 


14 


100 


2 


0.005 



THERA, to pp -> CD + X at the Tevatron (Run II), and to pp — > CD + X at the LHC, 
respectively. In each figure, parts (a) and (b) give the pr distributions da/dp T and the 
yc distributions da /dye, respectively. In each part, there are four frames, which refer to 
D = 7, Z, W, H, respectively. In each frame, we separately consider C = J/ip,XcJ, both 
in the CSM and in NRQCD. It is summed over C = Xco, Xa, Xc2 and D = W + ,W~. In 
the case of D = 7, the yc distributions are evaluated imposing the cut pr > 1 GeV in 
order to exclude the infrared and collinear singularities mentioned in Sec. 0. Figure |3](b) 
refers to the laboratory frame, ?/Jf b being related to its counterpart yc in the cm. frame 
byEq. 

We start the discussion of the figures with a few general observations. 

1. In all considered types of experiments, the associated production of Xcj+1, J / ip+W, 
Xcj + W, and J /if) + H is forbidden in the CSM to the order considered. 

2. The cross section ratio of J /ip + W and Xcj + W associated production, which pro- 



ceeds through processes (|T|) and (|T7|), is always (O j ^ 3 S[ 8) ) / £j=o lO XcJ 3 S| 



7(8) 



0.21 



With increasing value of mcD, the CO process (15) with q = 3 Si generally gains 
relative importance, since its cross section involves a gluon propagator with small 
virtuality, q 2 = m c , and is, therefore, enhanced by powers of m CD /m c relative to 
those of the other contributing processes. In the fragmentation picture f35| , this 
cross section would be evaluated by convoluting the cross section of qq — > gN with 



the g — >• cc 3 S[ fragmentation function [|3~6[1 . Processes fll6|) and (^7\ ) also benefit 
from this type of enhancement. 



The processes with D = 7 generally have much larger cross sections than those 
with D = Z, W, H because the available phase space is considerably ampler and the 
infrared and collinear singularities at px = 0, albeit eliminated by a minimum-p^ 
cut, still feed into the finite parts of the cross sections. 
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5. As mentioned in the Introduction, the cross sections of the processes with D = H 
are suppressed relative to those with D = Z, W by the smallness of the charm 
Yukawa coupling. 

6. As is evident from Eqs. (|TT|) and (|12|) , the direct-photoproduction channels in e + e~ 
and pe ± collisions correspond to pure CS and CO processes, respectively. 

7. If mcD 3> m Ci which is the case if pt 3> mp or D = Z, W, H, the resolved channels 
are generally suppressed against the direct ones because, according to Eq. (f|), the 
values of x a ,b are then close to unity, where the photon PDFs take small values. 

8. Direct photons participate in the hard scattering with their full momenta, while re- 
solved ones pass on only a fraction of theirs. Thus, in photoproduction at THERA, 
the direct and resolved cross sections are peaked in the backward and forward di- 
rections, respectively. 

We now discuss the e + e~ mode of TESLA [see Figs. 0(a) and (b)]. Here, J/ip + 7 
associated production dominantly proceeds through the direct CS process (|TT|) with q = 
3 S\ and iV = 7. Consequently, in the first frames of Figs. |l](a) and (b), the dotted and 
solid lines are superjacent. Xcj + 7 associated production only proceeds through singly- 
resolved or doubly-resolved CO processes and is accordingly suppressed relative to J /ip+1 
associated production. For the reasons exposed in the preceding paragraph, the reactions 
with D = Z chiefly proceed through the direct CS processes (|TT| ) with q = 3 Si, 3 Pj and 
N = Z. The contributing CO processes are singly or doubly resolved and, therefore, 
strongly suppressed. The reactions with D = W are mediated by the doubly-resolved CO 
processes ([Tj]) and (|17D and are accordingly suppressed. Xcj + H associated production 
dominantly proceeds through direct CS processes, while the contributing CO processes are 
doubly resolved and, therefore, heavily suppressed. Consequently, in the fourth frames 
of Figs. |l|(a) and (b), the short-dashed and medium-dashed lines are superjacent. On 
the other hand, J/ip + H associated production only proceeds through singly-resolved or 
doubly- resolved CO processes and is accordingly suppressed relative to Xcj + H associated 
production. 

We now turn to the 77 mode of TESLA [see Figs. ||(a) and (b)]. In contrast to 
the energy spectra of bremsstrahlung and beamstrahlung, which are strongly peaked 
at x = 0, the one of the back-scattered laser photons is evenly spread in the lower x 
range and exhibits a maximum at x = x max . According to Eq. (Q), large values of pt 
or mi) entail large values of x±. This explains why, in the 77 mode of TESLA, the pt 
spectra are less steep and the cross sections for D = Z,W, H are larger than in the e + e~ 
mode. Furthermore, the influence of the singly and doubly resolved channels is generally 
increased because the average photon energy is larger. Here, the CSM prediction for 
J/ijj + 7 associated production is dominated by the doubly-resolved process (|T3| ) with 
<^ = 3 5i and iV = 7 in the lower pt range, for pt ^ 5 GeV, while the direct process ( |TT|) 
with q = 3 Si and iV = 7 preponderates for larger values of px- Beholding the first frame of 
Fig. ||(b), we observe that the yc distribution of Xcj + 7 associated production is peaked 
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in the very forward and backward directions, close to the kinematic boundaries. This 
may be traced to the single-resolved process fll2|) with q = 3 Si and N = 7. This partonic 
subprocess also generates the pronounced shoulders in the yc distribution of J/ip + 7 
associated production in NRQCD, which is shown in the same figure. In the case of the 
e + e~ mode, it is suppressed by the softness of the effective photon-energy spectrum. The 
increased influence of the singly-resolved and doubly-resolved CO processes is reflected 
in the first frame of Fig. §(a) by the dispartment of the dotted and solid lines. Similar 
separations of CSM and NRQCD results are visible in the second frames of Figs. 0(a) and 
(b) for D = Z. 

We now move on to THERA [see Figs. |3|(a) and (b)]. In the CSM, <//■?/> + 7 associated 
production now only happens in the resolved channel, through process (|13|) with q = 3 Si 
and iV = 7, while, in NRQCD, it also takes place in the direct one, through process 
(|T2|) with q = 3 S\ and N = 7. This explains why the dotted line in the first frame of 
Fig. |3|(b) is peaked in the forward direction, while the solid one also exhibits a shoulder in 
the backward direction. The fact that this shoulder is not prominent may be understood 



by observing that (<D J ^ [ 3 Sr } \ ) / (p J ^ [ 3 Sj j J ) « 3.4 x 1(T 3 0. Detailed inspection 
reveals that, contrary to naive expectations, the direct contribution is at least one order 
of magnitude suppressed against the resolved one, even at large values of px- From the 
same figure, we also glean that the CSM result almost exhausts the resolved contribution 
in the lower px range, which is decisive for the yc distributions. In the intermediate px 
range, 6 <^px &25 GeV, the latter is dominated by process (|H|) with q = ^o, 3 Pj and 
N = 7, while, for yet larger values of px, process (|T5|) with q = 3 S\ and N = 7 gets in the 
lead.Q In the case of XcJ+7 associated production, processes ( |T2"| ) and ( |T5D with q = 3 S\ and 
N = 7 compete with each other. The former dominates for px ^ 11 GeV, while the latter 
prevails in the complementary px range. Since the yc distributions reflect the situation in 
the small-pr range, this explains why the medium-dashed line in the first frame of Fig. |3](b) 
is peaked in the backward direction. Also in the case of J/ip + Z associated production, 
the direct and resolved channels compete with each other. The direct one is dominated 
by process with q = ^o, 3 Pj and N — Z, while the resolved one is dominated by 
process ([T^) with q = 3 Si and N = Z. Here, the cross over occurs at p T ~ 28 GeV. By 
contrast, Xcj + Z associated production is always overwhelmingly dominated by the very 
same resolved process. This explains why the solid and medium-dashed lines in the second 
frame of Fig. 0(b) are peaked in the backward and forward directions, respectively. In the 
large-pr limit, the J/i/j + Z to Xcj + Z cross section ratio approaches the same value as in 
the case of D = W discussed under Item 2 in the above enumeration, processes ([T2]) The 
CS processes (|13|) with q = 3 Si, 3 Pj and N = Z are dramatically suppressed because they 
are resolved and not fragmentation prone. J ftp + H associated production dominantly 
proceeds through direct CO processes. On the other hand, Xcj + H associated production 
is only possible through resolved processes and is, therefore, suppressed relative to J / tp+H 
associated production. There are no fragmentation-prone processes in the case D = H 
because the Higgs boson always couples to the charm-quark line. 



3 This observation does not support the assertion made in Ref. [[l5| that light-quark-initiated processes 
are strongly suppressed and can be safely neglected. 
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We now proceed to the Tevatron [see Figs. f|(a) and (b)]. For pt^5 GeV, J/ip + 7 
associated production dominantly proceeds through the CS process ( |i~3"D with ^ = 3 Si and 



N = 7, due to the large value of (O j ^ 



3 Si ). For larger values of pt, the CO processes 



take over. Specifically, the latter are dominated by process (0) with ^ = ^o, 3 Pj and 
N = 7 for pt ^ 46 GeV, while process (|15|) with q = 3 Sx and iV = 7 prevails for larger 
values of pt- On the other hand, Xcj + 7 associated production is always dominated 
by the latter. Similarly, in the case of D = Z, process ([15]) with q = 3 Si and N = Z 
is always dominant. By the same token, the J/ip + Z to Xcj + Z cross section ratio is 
approximately the same as in the case of D = W discussed above. On the other hand, 
the CS processes ( |13"D with ^ = 3 Si, 3 Pj and N = Z are not fragmentation prone and are, 
therefore, vigorously suppressed. 

Finally, we arrive at the LHC [see Figs. ||(a) and (b)]. Here, the pt and yc distri- 
butions generally exhibit very similar shapes as in the case of the Tevatron, while their 
normalizations are increased by approximately the same amount as the value of v^S" is. 

The experimental observation of processes that are forbidden or exceedingly suppressed 
in the CSM with cross sections that are compatible with the NRQCD predictions would 
provide striking evidence for the NRQCD factorization hypothesis and firmly establish 
the existence of CO processes in nature. Furthermore, precise measurements of these 
cross sections would lead to useful constraints on the appearing CO MEs. On the other 
hand, the experimental study of processes that are dominated by CS channels would allow 
for independent determinations of the CS MEs, which can be compared with the results 
obtained using traditional methods, so as to lead to valuable consistency checks. Our 
analysis allows us to identify processes of all these categories. As mentioned above, the 
associated production of Xcj + 7, J/ip + W, Xcj + W, and J/ip + H exclusively proceeds 
through CO processes. At THERA, the Tevatron, and the LHC, the associated production 
of J/i[) + r y and Xcj + H with px 3> mc and of J /ip + Z and Xcj + Z with arbitrary values 
of pt is greatly dominated by CO processes. On the other hand, examples where CO 
processes play a negligible role include the associated production of J/tj) + 7, J /if) + Z, 
Xcj + Z, and Xcj + H in the e + e~ mode of TESLA, of Xcj + H in the 77 mode of TESLA, 
and of J/ip + 7 and Xcj + H with small values of a t THERA, the Tevatron, and the 
LHC. 

We conclude this section by assessing the observability of the various processes in the 
various experiments with the aid of Table |I[ The processes with D = 7 will abundantly 
take place in all considered experiments. The processes with D = Z,W will produce 
considerable yields at the hadron colliders, namely, several hundred (ten thousand) events 
per year at the Tevatron (LHC), while they significantly fall short of the one-event-per- 
year mark at TESLA and THERA. As expected, the processes with D = H are predicted 
to be far too rare to be observable in any of the considered experiments. In turn, their 
observation would hint at new physics beyond the SM. 
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4 Conclusions 



We studied the associated production of heavy quarkonia C with electroweak bosons D 
in photon-photon, photon-hadron, and hadron-hadron collisions to LO in the NRQCD 
factorization formalism. We considered all experimentally established heavy quarkonia, 
with 2S+1 Lj = 1 5' , 3 5i, 1 .P 1 , 3 .P i7 , and all electroweak bosons of the SM, D = -f,Z,W,H. 
We listed all contributing partonic cross sections, except for those with D = H, which 
will be irrelevant for phenomenology in the foreseeable future. We presented numerical 
results for any combination of C — J/ip, Xcj and D = 7, Z, W, H appropriate for TESLA 
in the e + e~ and 77 modes, THERA, Run II of the Tevatron, and the LHC. 

At TESLA and THERA, only the processes with D = 7 will have observable cross 
sections, while, at the Tevatron and the LHC, this is also the case for D = Z, W. Ob- 
servation of Xcj + 7, J/i> + W, and Xcj + W associated production, which are pure CO 
processes, would give strong support to the NRQCD factorization hypothesis and, if mea- 
sured with sufficient precision, allow for independent determinations of (o J l'^ 3 s[ 8 ^ ^ and 

(O Xc0 3 Sj 8 ' ) "j. A similar statement applies to J/ip + Z and Xcj + Z associated production 
at the Tevatron and the LHC, which are overwhelmingly dominated by CO processes. 
On the other hand, CO processes play a subordinate role for J/ip + 7 associated pro- 
duction in the e + e~ mode of TESLA for arbitrary values of px and at THERA, the 
Tevatron, and the LHC for small values of pt- This offers the opportunity to extract 
a new value of (o J ^ 3 S^ ^ and thus to allow for a consistency check. Alternatively, 

assuming (o J ^ 3 Si ^ to be sufficiently well known, one may get a better handle on the 
gluon PDFs of the photon and proton by fitting THERA, Tevatron, and LHC data of 
7/^ + 7 associated production in the lower pt range. 
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A Partonic cross sections 

In this appendix, we list the differential cross sections da/dt for processes (|ll"|)-(|i~7|) with 
N = Z to be inserted in Eq. (|l|) and its counterparts for lepton-hadron and hadron-hadron 
scattering. The results for N = 7 are recovered as explained in Sec. 0. Our results read 
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- s(1741t 6 + 7835t 5 M + 10749t 4 M 2 + 8966t 3 M 3 + 10749tV + 7835to 5 + 1741m 6 ) 

- (t + u)(99t 6 + 343t 5 M - 679tV - 2118tV - 679t 2 M 4 + 343tM 5 + 99m 6 )] 
+ m|[584s 8 + 2956s 7 (t + u) + 4s 6 (1259t 2 + 2842£m + 1259m 2 ) 

+ 2s 5 (t + w)(237t 2 + 2626tw + 237m 2 ) 

- s 4 (9273t 4 + 30930t 3 M + 40826* V + 30930tw 3 + 9273m 4 ) 

- 2s 3 (t + M)(6289t 4 + 22639t 3 M + 28320tV + 22639£m 3 + 6289m 4 ) 

- s 2 (7039t 6 + 40430t 5 M + 86805* V + 105964* V 
+ 86805t 2 M 4 + 40430tM 5 + 7039m 6 ) 

- 2s(t + m)(811* 6 + 4992t 5 M + 9574tV + 9982* V + 9574tV + 4992*m 5 + 811m 6 ) 

- 2(39* 8 + 404t 7 M + 1209tV + 1488* 5 m 3 + 1280tV 
+ 1488t 3 M 5 + 1209t 2 M 6 + 404*m 7 + 39m 8 )] 

- 2m|[56s 9 + 386s 8 (* + u) + 2s 7 (539* 2 + 1164*m + 539m 2 ) 
+ 2s 6 (t + m)(595* 2 + 1466*m + 595m 2 ) 

- 2s 5 (201* 4 + 599t 3 M + 588* V + 599*m 3 + 201m 4 ) 

- 2s 4 (* + m)(1189* 4 + 4965* 3 m + 6653*V + 4965*m 3 + 1189m 4 ) 

- s 3 (2455* 6 + 15774* 5 m + 37813* 4 m 2 + 48736* 3 m 3 + 37813* 2 m 4 + 15774*m 5 + 2455m 6 ) 

- s 2 (* + m)(1127* 6 + 7744* 5 m + 18069* V + 22220*V 
+ 18069* 2 m 4 + 7744*m 5 + 1127m 6 ) 

- s(213* 8 + 2150* 7 m + 7654*V + 14044* 5 m 3 + 16614*V 
+ 14044*V + 7654*V + 2150*m 7 + 213m 8 ) 

- (* + u)(7t 8 + 118* 7 m + 528* V + 978* 5 m 3 + 1098*V 
+ 978*V + 528* V + 118*m 7 + 7m 8 )] 

+ m A z [8s 10 + 76s 9 (* + u) + 8s 8 (44* 2 + 93*m + 44m 2 ) 

+ 6s 7 (* + m)(127* 2 + 278*m + 127m 2 ) 

+ s 6 (597* 4 + 2324* 3 m + 3574* V + 2324*m 3 + 597m 4 ) 

- 2s 5 (* + m)(233* 4 + 1377* 3 m + 1950* V + 1377*m 3 + 233m 4 ) 

- s 4 (1327* 6 + 9644* 5 m + 25033*V + 33280* 3 m 3 + 25033*V + 9644*m 5 + 1327m 6 ) 

- 2s 3 (* + m)(547* 6 + 4136* 5 m + 10580*V + 13738*V 
+ 10580* 2 m 4 + 4136*m 5 + 547m 6 ) 

- s 2 (413* 8 + 4348* 7 m + 16425* V + 32576* V + 40108*V 
+ 32576*V + 16425*V + 4348*m 7 + 413m 8 ) 

- 2s(t + m)(31* 8 + 409* 7 m + 1674*V + 3332*V + 4016* 4 m 4 
+ 3332*V + 1674* 2 m 6 + 409*m 7 + 31m 8 ) 

- (* + u) 2 {f + 34* 7 m + 225*V + 524*V + 628*V 
+ 524* 3 m 5 + 225*V + 34*m 7 + m 8 )] 
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- 2m|[6s 9 (t + uf + 30s 8 (t + uf + s 7 {t + m) 2 (5K 2 + 83tu + 51m 2 ) 
+ s 6 (t + m)(13£ 4 - 39t 3 M - 98tV - 39£m 3 + 13m 4 ) 

- s 5 (70t 6 + 613t 5 M + 1782t 4 M 2 + 2470tV + 1782tV + 613tn 5 + 70m 6 ) 

- s 4 (t + m)(104£ 6 + 873t 5 M + 2434tV + 3298* V + 2434* V + 873*m 5 + 104m 6 ) 

- s 3 (65* 8 + 733* 7 m + 2918* 6 m 2 + 6080tV + 7652* V 
+ 6080* 3 m 5 + 2918* 2 m 6 + 733*m 7 + 65m 8 ) 

- s 2 (* + m)(19* 8 + 263* 7 m + 1082* V + 2252* 5 m 3 + 2802* V 
+ 2252*V + 1082*V + 263*m 7 + 19m 8 ) 

- s{t + m) 2 (2* 8 + 43* 7 m + 200* V + 440*V + 540* V 
+ 440* 3 m 5 + 200*V + 43*m 7 + 2m 8 ) 

- tu(t + uf(t 2 + tu + m 2 )(* 4 + 11* 3 m + 24*V + 11*m 3 + m 4 )] 

- (t + u)(s + t + m) 2 [s 6 (* + u)(t 2 + 4*m + u 2 ) 
+ 2s 5 (2* 4 + m 3 M + 20* V + 11*m 3 + 2m 4 ) 

+ 2s 4 (* + m)(3* 4 + 17* 3 m + 32*V + 17*m 3 + 3m 4 ) 

+ 2s 3 (* + m) 2 (2* 4 + 15* 3 m + 24*V + 15*m 3 + 2m 4 ) 

+ s 2 (* + m)(* 6 + 18* 5 m + 48*V + 68* 3 m 3 + 48*V + 18*m 5 + m 6 ) 

+ 2s*m(* 2 + tu + m 2 )(2* 4 + 7* 3 m + 12*V + 7*m 3 + 2m 4 ) + *V(* + m) 5 ]}, (A.6) 



F 2 = 96mf - 64mf [lis + 15(* + u)} 

+ 8m|°[192s 2 + 713s(* + u) + 28(17* 2 + 35*m + 17m 2 )] 

+ 4m^ 8 [40s 3 - 2689s 2 (* + u) - 2s(2359* 2 + 4960*m + 2359m 2 ) 

- 12(* + m)(173* 2 + 384*m + 173m 2 )] 

- 2m^ 6 [3020s 4 - 610s 3 (* + u) - 2s 2 (7615* 2 + 16398*m + 7615m 2 ) 

- 3s(* + m)(5775* 2 + 13402*m + 5775m 2 ) 

- 5647* 4 - 25540* 3 m - 39866*V - 25540*m 3 - 5647m 4 ] 

+ 4mi 4 [2836s 5 + 5845s 4 (* + u) - 2s 3 (730* 2 + 1927*m + 730m 2 ) 

- 2s 2 (* + m)(5965* 2 + 14477*m + 5965m 2 ) 

- s(9956* 4 + 46445* 3 m + 73390*V + 46445*m 3 + 9956m 4 ) 

- 2(* + m)(1269* 4 + 6380* 3 m + 10334*V + 6380*m 3 + 1269m 4 )] 

- m% [10600s 6 + 36980s 5 (* + u) + 4s 4 (9559* 2 + 19444*m + 9559m 2 ) 

- 2s 3 (* + m)(5559* 2 + 16706*m + 5559m 2 ) 

- s 2 (47103* 4 + 224004* 3 m + 357674* V + 224004*m 3 + 47103m 4 ) 

- 2s(t + m)(15175* 4 + 78107* 3 m + 129108* V + 78107*m 3 + 15175m 4 ) 

- 6165t 6 - 48348* 5 m - 139987*V - 195640*V 

- 139987*V - 48348*m 5 - 6165m 6 ] 
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+ 2m^°[2880s 7 + 13554s 6 (t + u) + 4s 5 (6177t 2 + 12874tM + 6177m 2 ) 
+ 2s 4 (t + w)(8249t 2 + 18092tM + 8249m 2 ) 

- s 3 (6683t 4 + 31527t 3 M + 51696tV + 31527tM 3 + 6683m 4 ) 

- s 2 (t + w)(15723t 4 + 78965t 3 M + 132576* V + 78965*m 3 + 15723m 4 ) 

- 3s(2621t 6 + 2041 lt 5 w + 59773tV + 84038* 3 m 3 
+ 59773tV + 20411tw 5 + 2621m 6 ) 

-(t + M)(1257t 6 + 11553t 5 M + 36251tV + 51926t 3 M 3 
+ 3625UV + 11553tM 5 + 1257m 6 )] 

- m| [1864s 8 + 10820s 7 (t + m) + 4s 6 (6847t 2 + 14418tM + 6847m 2 ) 
+ 2s 5 (t + w)(17033t 2 + 39274tw + 17033m 2 ) 

+ s 4 (14141t 4 + 71138t 3 M + 111826tV + 71138tM 3 + 14141m 4 ) 

- 2s 3 (t + M)(5695t 4 + 21723t 3 M + 37444tV + 21723£m 3 + 5695m 4 ) 

- s 2 (14753t 6 + 104550t 5 M + 300227* V + 421692t 3 M 3 
+ 300227t 2 M 4 + 104550tM 5 + 14753m 6 ) 

- 2s(t + m)(277U 6 + 23944t 5 M + 75106t 4 M 2 + 108454t 3 M 3 
+ 75106* V + 23944tM 5 + 2771m 6 ) 

- 2(t + m) 2 (331* 6 + 3842t 5 M + 13494* 4 m 2 + 19926* 3 m 3 
+ 13494tV + 3842tM 5 + 331m 6 )] 

+ 2m 6 [168s 9 + 1158s 8 (t + u) + 12s 7 (316t 2 + 671tM + 316m 2 ) 
+ 4s 6 (t + M)(1693t 2 + 3990tM + 1693m 2 ) 
+ s 5 (5796t 4 + 29542t 3 M + 47188tV + 29542tM 3 + 5796m 4 ) 
+ 2s 4 (t + M)(162t 4 + 4996t 3 M + 8301tV + 4996tM 3 + 162m 4 ) 

- s 3 (3347t 6 + 16010t 5 M + 39319tV + 53548t 3 M 3 
+ 39319tV + 16010tM 5 + 3347m 6 ) 

- s 2 (t + M)(2439t 6 + 16826t 5 M + 49619t 4 M 2 + 71420£ 3 m 3 
+ 49619tV + 16826tM 5 + 2439m 6 ) 

- s(645t 8 + 7666t 7 M + 35040t 6 M 2 + 82140tV + 108250t 4 M 4 
+ 82140tV + 35040* V + 7666tM 7 + 645m 8 ) 

- (t + M) 3 (51t 6 + 852t 5 M + 3539tV + 5396tV + 3539tV + 852£m 5 + 51m 6 )] 

- m|[24s 10 + 212s 9 (t + u) + 4s 8 (24K 2 + 512£m + 241m 2 ) 
+ 2s 7 (t + M)(1179t 2 + 2822tM + 1179m 2 ) 

+ s 6 (2855t 4 + 15364t 3 M + 24930tV + 15364tM 3 + 2855m 4 ) 
+ 2s 5 (t + M)(367t 4 + 6085t 3 M + 10618tV + 6085£m 3 + 367m 4 ) 

- s 4 (2117t 6 + 1876t 5 M - 6829tV - 13152t 3 M 3 - 6829tV + 1876£m 5 + 2117m 6 ) 

- 2s 3 (t + m)(124U 6 + 4958t 5 M + 11056tV + 15114t 3 M 3 
+ 11056t 2 M 4 + 4958tM 5 + 1241m 6 ) 
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- s 2 (1079t 8 + 9328t 7 M + 36475tV + 80988t 5 M 3 + 105492t 4 M 4 
+ 80988tV + 36475t 2 M 6 + 9328tM 7 + 1079m 8 ) 

- 2s(t + m)(9K 8 + 1293t 7 M + 6382tV + 15744t 5 M 3 + 21112tV 
+ 15744tV + 6382t 2 M 6 + 1293tM 7 + 91m 8 ) 

- (t + u) 4 (7t 6 + 228£ 5 m + 1252tV + 1952t 3 M 3 + 1252tV + 228tM 5 + 7m 6 )] 
+ 2m|[16s 9 (t + uf + 2s 8 (t + w)(35t 2 + S2tu + 35m 2 ) 

+ s 7 (95t 4 + 627t 3 M + 1056tV + 627£m 3 + 95m 4 ) 

- s 6 (t + u)(57t 4 - 721t 3 M - 1422iV - 721£m 3 + 57m 4 ) 

- s 5 (332t 6 + 227t 5 M - 1538* V - 2874tV - 1538t 2 M 4 + 227£m 5 + 332m 6 ) 

- s 4 (t + w)(424t 6 + 911t 5 M + 608tV + 298tV + 608tV + 911tw 5 + 424m 

- s 3 (253t 8 + 1415t 7 M + 3802tV + 6836* V + 8376* V 



+ 6836* 3 m 5 + 3802tV + 1415*m 7 + 253m 8 ) 
- s 2 (t + m)(69* 8 + 61K 7 m + 2282*V + 5134* V + 6750*V 
+ 5134*V + 2282* 2 m 6 + 611*m 7 + 69m 8 ) 
-s(* + m) 2 (6* 8 + 129* 7 m 



682*V + 1796* 5 m 3 + 2452* 4 m 



5„.3 



+ 1796*V + 682* 2 m 6 + 129*m 7 + 6m 



- tu(t + m) 5 (7* 4 + 68* 3 m + 105* V + 68*m 3 + 7m 4 )] 
+ (t + u)(s + t + m) 2 [s 6 (* + u)(7t 2 - Atu + 7m 2 ) 
+ 2s 5 (14* 4 + 9* 3 m - 4*V + 9*m 3 + 14m 4 ) 
+ 2s 4 (* + m)(21* 4 + 3* 3 m - 4*V + 3*m 3 + 21m 4 ) 
+ 2s 3 (* + m) 2 (14* 4 + 5* 3 m + 10*V + 5*m 3 + 14m 4 ) 
+ s 2 (t + m)(7* 6 + 38* 5 m + 64* 4 m 2 + 116* 3 m 3 + 64* V + 38*m 5 + 7m 6 ) 
+ 2s*m(6* 6 + 25* 5 m 



69*V 



94*V + 69*V + 25*m 5 + 6m 6 ) 



x2..4 



+ *V(* + m) 3 (7* 2 + lOtu + 7m 2 )], 



(A.7) 



dt 



-(i) 



? = 1 S'o, 3 S , i, 1 Pi, (A.8) 



3p(8) 



da 



J'=0 



(gg QQ[n]Z) = ^ 4 % (77 - QQ[n]z) , n = 3 S[", 'P« 'Pj 1 ', 



dt 



(A.10) 
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KaW 

12Mm|s 3 (m| - s - t) 2 (m 2 z -s- u) 2 {2m 2 z - t - u) 2 

x {9wjm|[9m^ - 3mf(8s + 9(t + u)) 

+ m|(24s 2 + 56s(t + u) + 3(llt 2 + 23tu + 11m 2 )) 

- m|(10s 3 + 42s 2 (t + u) + 2s(25t 2 + 54tu + 25m 2 ) + 3(f + u)(7t 2 + 16fu + 7m 2 )) 
+ m z (s 4 + 12s 3 (t + u) + s 2 (25t 2 + 58£m + 25m 2 ) + As(t + u)(5t 2 + Utu + 5m 2 ) 

+ 7t 4 + 35t 3 M + 51* V + 35tM 3 + 7m 4 ) 

- m|(s 4 (t + m) + 2s 3 (2t 2 + 5tM + 2m 2 ) + 5s 2 (t + u){t 2 + Atu + m 2 ) 
+ s(t + w) 2 (3t 2 + 16tu + 3m 2 ) + (t + u)(t 2 + tu + u 2 ){t 2 + 7tu + m 2 )) 
+ tu(s 4 + 2s 3 (t + u) + 3s 2 (t + m) 2 + 2s(t + m) 3 + (t 2 + tu + u 2 ) 2 )] 

+ ha 2 Q s(2m 2 z - s-t- u) 2 [m% - 2m%(t + u) - m%(s 2 - t 2 - Atu - u 2 ) 

-2m 2 z tu(t + u) +t 2 u 2 ]}, (A.11) 



QQ 



lo(8) 
D 



36Mm|s 3 (m| - s - t) 2 (m 2 z - s - u) 2 {2m 2 z - t - u) 2 
x {hv 2 Q m 2 z s[m l £ - 4m|(3s + t + u) + m|(22s 2 + 26s(t + «) + 5t 2 + 12fu + 5m 2 ) 

- 2m|(5s 3 + 14s 2 (t + m) + s(8t 2 + 23£m + 8m 2 ) + (i + M)(t 2 + 5£m + m 2 )) 

- m|(s 4 - 4s 3 (t + u) - s 2 (9t 2 + 26tM + 9m 2 ) - 2s(t + M)(t 2 + lOtu + m 2 ) 

-tM(4t 2 + 9tM + 4M 2 )) 

- 2(s 3 (t 2 + tu + u 2 ) + s 2 {t + m) 3 + stM(t 2 + 3tu + u 2 ) + t 2 u 2 {t + m))] 

- 9aQ[9m^ 2 (s + t + u) - m^°(18s 2 + 49s(t + u) + 27(t + m) 2 ) 

+ m|(7s 3 + 57s 2 (t + u) + s(83t 2 + 169£m + 83m 2 ) + 3(t + M)(llt 2 + 23tM + 11m 2 )) 
+ m|(2s 4 - 13s 3 (t + u) - 2s 2 (31t 2 + 63£m + 31m 2 ) - 5s(t + M)(13t 2 + 28£m + 13m 2 ) 

- 3(t + M) 2 (7t 2 + 16fu + 7m 2 )) 

- m 4 z {2s 4 {t + u) - s 3 (8t 2 + 19tM + 8m 2 ) - 5s 2 (t + u)(3f + 2m) (2t + 3m) 

- 5s(5t 4 + 23t 3 M + 35t 2 M 2 + 23tM 3 + 5m 4 ) 

- (t + u)(7t 4 + 35t 3 M + 51t 2 M 2 + 35tM 3 + 7m 4 )) 

+ m|(s 4 (t 2 + m 2 ) - 2s 3 (t + M)(t 2 + 3tM + m 2 ) - s 2 (t + M) 2 (6t 2 + Yltu + 6m 2 ) 

- s(t + M)(4t 4 + 26t 3 M + 37t 2 M 2 + 26tM 3 + 4m 4 ) 

- (t + uf(t 4 + 8t 3 M + 9t 2 M 2 + 8tM 3 + M 4 )) 

+ tu(s + t + u)(s(t + u)+t 2 + tu + u 2 ) 2 ]}, (A.12) 



7/7 



99 - QQ 6 S± 



3c(8) 



da / 
dt(" 



QQ [^J Z) 



-7ia 2 g 2 



9M 3 m|s 3 (m| - s - tf{m 2 z -s- m) 3 (2tm| -t-u) 
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x {10^m|s[48m^ 8 - 16m^ 6 (lls + 15(t + «)) 

+ 2mi 4 (141s 2 + 374s(t + u) + 4(65t 2 + 134tM + 65m 2 )) 

- 4m^ 2 (65s 3 + 255s 2 (t + u) + s(341t 2 + 696£m + 341m 2 ) 
+ 32(t + u)(5t 2 + l\tu + 5m 2 )) 

+ m^°(138s 4 + 798s 3 (t + u) + 10s 2 (158t 2 + 315tM + 158m 2 ) 

+ s(t + w)(1399t 2 + 2958£m + 1399m 2 ) 

+ 491t 4 + 2188t 3 M + 3378fV + 2188tM 3 + 491m 4 ) 

- m|(24s 5 + 356s 4 (t + u) + 20s 3 (53t 2 + lOltu + 53m 2 ) 

+ 10s 2 (t + M)(138t 2 + 265tM + 138m 2 ) + s(t + m) 2 (887£ 2 + 1956tw + 887m 2 ) 
+ (t + M)(241t 4 + 1188t 3 M + 1846tV + 1188tM 3 + 241m 4 )) 

- m|(18s 6 - 60s 5 (t + u) - s 4 (417t 2 + 734tM + 417m 2 ) 

- s 3 (t + M)(801t 2 + 1298tM + 801m 2 ) 

- s 2 (749t 4 + 2772t 3 M + 4068tV + 2772tM 3 + 749m 4 ) 

- s(t + M)(359t 4 + 1534t 3 M + 2362tV + 1534£m 3 + 359m 4 ) 

- 2(t + M) 2 (37t 4 + 217t 3 M + 332tV + 217tM 3 + 37m 4 )) 
+ m 4 z (12s 7 + 16s 6 (t + u) - s 5 (83t 2 + \22tu + 83m 2 ) 

- s 4 (t + M)(275t 2 + 356tM + 275m 2 ) 

- s 3 (367t 4 + 1153t 3 M + 1580tV + 1153tM 3 + 367m 4 ) 

- s 2 (t + m)(258£ 4 + 843t 3 M + 1214tV + 843£m 3 + 258m 4 ) 

- 4s(23t 6 + 143t 5 M + 371fV + 500t 3 M 3 + 37UV + 143£m 5 + 23m 6 ) 
-(t + m) 3 (13£ 4 + 104t 3 M + 150tV + 104tM 3 + 13m 4 )) 

- m|(2s 8 + 6s 7 {t + u) - s G {7t 2 + 4tu + 7m 2 ) - 52s 5 (t + M)(t 2 + tu + m 2 ) 

- s 4 (97t 4 + 269t 3 M + 346tV + 269tM 3 + 97m 4 ) 

- 2s 3 (t + M)(48t 4 + 122t 3 M + 155tV + 122tM 3 + 48m 4 ) 

- s 2 (53t 6 + 244t 5 M + 531tV + 674t 3 M 3 + 531tV + 244£m 5 + 53m 6 ) 

- s(t + M)(14t 6 + 82t 5 M + 217* V + 290t 3 M 3 + 217tV + 82tM 5 + 14m 6 ) 

- (t + M) 4 (t 4 + 15t 3 M + 19tV + 15tM 3 + m 4 )) 

- (t + u)(s 7 (t + u) + s 6 (5t 2 + 6tu + 5m 2 ) + s 5 (t + M)(llt 2 + Stu + 11m 2 ) 
+ 2s 4 (7t 4 + 16t 3 M + 19tV + 16tM 3 + 7m 4 ) 

+ s 3 (t + M)(m 4 + 21t 3 M + 25tV + 21tM 3 + 11m 4 ) 

+ s 2 (5t 6 + 19t 5 M + 38iV + 46t 3 M 3 + 38tV + 19tM 5 + 5m 6 ) 

+ s(t + M)(t 6 + 5t 5 M + 14t 4 M 2 + 18t 3 M 3 + 14* V + 5tM 5 + m 6 ) 

+ tM(t + M) 4 (t 2 +tM + M 2 ))] 

+ 9a 2 Q [36mf(s - t - u) - Amf(17s 2 + Ss(t + u) - 45(t + m) 2 ) 
+ m^ 6 (28s 3 + 12s 2 (t + u) - s(257t 2 + 466tM + 257m 2 ) 
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- 3(t + u)(131t 2 + 278tu + 131m 2 )) 

- 4mi 4 (2s 4 - 40s 3 (t + u) - 6s 2 (22t 2 + Altu + 22m 2 ) 

- 4s(f + m)(47£ 2 + 88tu + 47m 2 ) - 3(t + w) 2 (41t 2 + 98tu + 41m 2 )) 
+ m l i(28s b - 88s 4 (t + u) - s 3 (695t 2 + 1334£m + 695m 2 ) 

- s 2 (t + w)(1201t 2 + 2218tw + 1201m 2 ) 

- 2s(485t 4 + 1938t 3 M + 2910tV + 1938tM 3 + 485m 4 ) 

- 4(t + M)(97t 4 + 473t 3 w + 750tV + 473£m 3 + 97m 4 )) 

- m^°(20s 6 + 44s 5 (t + u) - 26s 4 (llt 2 + 20tM + 11m 2 ) 

- 8s 3 (t + M)(132t 2 + 241tM + 132m 2 ) 

- s 2 (1277t 4 + 4942t 3 M + 7386tV + 4942£m 3 + 1277m 4 ) 

- 2s(t + M)(363t 4 + 1544t 3 M + 2346* V + 1544tM 3 + 363m 4 ) 

- 6(t + M) 2 (33t 4 + 190t 3 M + 310tV + 190tM 3 + 33m 4 )) 
+ m|(4s 7 + 32s 6 (t + u) + s 5 (7t 2 + 30tM + 7m 2 ) 

- s 4 (t + M)(351t 2 + 566tM + 351m 2 ) 

- 2s 3 (421t 4 + 1577t 3 M + 2354tV + 1577£m 3 + 421m 4 ) 

- s 2 (t + M)(779t 4 + 3112t 3 M + 4734tV + 3112£m 3 + 779m 4 ) 

- 2s(168t 6 + 1138t 5 M + 2987t 4 M 2 + 4042t 3 M 3 + 2987t 2 M 4 + 1138tM 5 + 168m 6 ) 

- 2(t + M)(32t 6 + 303t 5 M + 915fV + 1280t 3 M 3 + 915tV + 303tM 5 + 32m 6 )) 

- m| (4s 7 (t + m) + 2s 6 (9t 2 + 16tu + 9m 2 ) - 4s 5 (t + m)(5£ 2 + 6tM + 5m 2 ) 

- 2s 4 (107t 4 + 374t 3 M + 550tV + 374£m 3 + 107m 4 ) 

- 3s 3 (t + M)(127t 4 + 478t 3 M + 742tV + 478£m 3 + 127m 4 ) 

- 2s 2 (141t 6 + 897t 5 M + 2309tV + 3134fV + 2309tV + 897tM 5 + 141m 6 ) 

- s(t + M)(95t 6 + 746t 5 M + 2017tV + 2764tV + 2017tV + 746£m 5 + 95m 6 ) 

- A{t + M) 2 (3t 6 + 40t 5 M + 134t 4 M 2 + 186fV + 134t 2 M 4 + 40tM 5 + 3m 6 )) 
+ m|(s 7 (t + m) 2 + s 6 (t + u)(5t 2 + 2tu + 5m 2 ) 

- 2s 5 (5t 4 + 24t 3 M + 36tV + 24tM 3 + 5m 4 ) 

- s 4 (t + M)(65t 4 + 240t 3 M + 374tV + 240£m 3 + 65m 4 ) 

- s 3 (94t 6 + 577t 5 M + 1489* V + 2060t 3 M 3 + 1489* V + 577*m 5 + 94m 6 ) 

- s 2 (t + m)(57* 6 + 424t 5 M + 1100t 4 M 2 + 1578t 3 M 3 + 1100*V + 424*m 5 + 57m 6 ) 

- s(15t 8 + 186t 7 M + 758tV + 1626fV + 2074tV 
+ 1626t 3 M 5 + 758t 2 M 6 + 186tM 7 + 15m 8 ) 

-(t + u)(t 8 + 28t 7 M + 160tV + 376* V + 490tV 
+ 376*V + 160t 2 M 6 + 28tM 7 + m 8 )) 

- m|(s 6 (t - u) 2 (t + m) 2 - s 5 (t + u)(t 4 + 12t 3 M + 18tV + 12tM 3 + m 4 ) 

- 2s 4 (4t 6 + 29t 5 M + 77* V + llOtV + 77* V + 29*m 5 + 4m 6 ) 

- 2s 3 (* + M)(5t 6 + 42t 5 M + 106t 4 M 2 + 170t 3 M 3 + 106t 2 M 4 + 42£m 5 + 5m 6 ) 
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- s 2 (5t 8 + 69t 7 M + 261tV + 577t 5 u 3 + 756t 4 M 4 + 577tV + 261t 2 M 6 + 69tu 7 + 5m 8 ) 

- s(t + w) 3 (t 6 + 19t 5 M + 53t 4 M 2 + 78fV + 53t 2 M 4 + 19tu 5 + m 6 ) 

- 2tu(t + w) 2 (t 6 + 9t 5 M + 21tV + 28tV + 21tV + 9tu 5 + m 6 )) 

- tu(t + u)(s 5 (t + u){t 2 + tu + u 2 ) + s 4 (4t 4 + 9t 3 M + 18tV + 9tu 3 + 4m 4 ) 
+ 2s 3 (t + u)(3t 4 + 5t 3 u + Ut 2 u 2 + 5£m 3 + 3m 4 ) 

+ 2s 2 (t + u) 2 (2t 4 + 3t 3 u + 9t 2 u 2 + 3tu 3 + 2m 4 ) 

+ S(t + M)(t 3 + 3t 2 M + 2tM 2 + M 3 )(t 3 + 2t 2 M + 3tM 2 + M 3 ) 

+ tM(t + M) 2 (t 2 +tM + M 2 ) 2 )]}, (A.13) 



™W 

3M 3 m|s 3 (m| - s - tf(m 2 z - s - u)\2m\ - t - m) 4 
x {9^Qm|[72m|° - 12m^ 8 (31s + 33(t + u)) 
+ m^ 6 (760s 2 + 1788s(t + u) + 975t 2 + 2046tM + 975m 2 ) 

- m^ 4 (804s 3 + 3144s 2 (t + u) + s(3785t 2 + 7854tM + 3785m 2 ) 
+ 108(t + M)(13t 2 + 30tM + 13m 2 )) 

+ m^ 2 (520s 4 + 2868s 3 (t + u) + s 2 (5661t 2 + 11578tM + 5661m 2 ) 

+ 4s(t + M)(1151t 2 + 2569tM + 1151m 2 ) 

+ 4(323t 4 + 1513t 3 M + 2376tV + 1513£m 3 + 323m 4 )) 

- m^°(236s 5 + 1644s 4 (t + u) + s 3 (4431t 2 + 8902tM + 4431m 2 ) 
+ 2s 2 (t + M)(2885t 2 + 6188tM + 2885m 2 ) 

+ s(3514t 4 + 15863t 3 M + 24534tV + 15863£m 3 + 3514m 4 ) 

+ 2(t + M)(389t 4 + 2086t 3 M + 3366tV + 2086£m 3 + 389m 4 )) 

+ m|(56s 6 + 628s 5 (t + u) + s 4 (2199t 2 + 4382£m + 2199m 2 ) 

+ 20s 3 (t + M)(190t 2 + 393tM + 190m 2 ) 

+ s 2 (3606t 4 + 15715t 3 M + 23922tV + 15715£m 3 + 3606m 4 ) 

+ s{t + M)(1720t 4 + 8657t 3 M + 13402t 2 M 2 + 8657£m 3 + 1720m 4 ) 

+ 2(150t 6 + 130H 5 m + 3823t 4 M 2 + 5328tV + 3823tV + 130Hm 5 + 150m 6 )) 

+ m|(4s 7 - 96s 6 (t + u) - 3s 5 (219t 2 + 442£m + 219m 2 ) 

- 4s 4 (t + M)(387t 2 + 800tM + 387m 2 ) 

- 4s 3 (479t 4 + 2045t 3 M + 3095tV + 2045£m 3 + 479m 4 ) 

- 2s 2 (t + M)(695t 4 + 3346t 3 M + 5022tV + 3346tM 3 + 695m 4 ) 

- s(526t 6 + 4237t 5 M + 11717t 4 M 2 + 15936t 3 M 3 + 11717tV + 4237tM 5 + 526m 6 ) 

- 4(t + M)(17t 6 + 198t 5 M + 636tV + 890t 3 M 3 + 636tV + 198tM 5 + 17m 6 )) 

- m|(4s 7 (t + u) - 69s 6 (t + m) 2 - 20s 5 (t + M)(17t 2 + 36tM + 17m 2 ) 

- 2s 4 (297t 4 + 1276t 3 M + 1954t 2 M 2 + 1276tM 3 + 297m 4 ) 



99 -> QQ 



3p 



(8) 
J 
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- 2s 3 (t + m)(276£ 4 + 1311t 3 M + 1996tV + Ulltu 3 + 276m 4 ) 

- s 2 (312t 6 + 2433t 5 M + 6563t 4 M 2 + 8836t 3 M 3 + 6563t 2 M 4 + 2433£m 5 + 312m 6 ) 

- s(t + u)(92t 6 + 947t 5 M + 2734£ 4 m 2 + 3626t 3 M 3 + 2734* V + 947tM 5 + 92m 6 ) 

- 7t 8 - 164t 7 M - 876tV - 2048tV - 2666tV 

- 2048tV - 876tV - lQAtu 7 - 7m 8 ) 

+ m 2 z (s 7 (t + uf - 2s 6 (t + u)(llt 2 + 28tu + 11m 2 ) 

- s 5 (80t 4 + 387t 3 M + 590tV + 387£m 3 + 80m 4 ) 

- 4s 4 (t + w)(26t 4 + 143t 3 M + 215tV + 143£m 3 + 26m 4 ) 

- s 3 (72t 6 + 635t 5 w + 1741tV + 2396tV + 1741tV + 635£m 5 + 72m 6 ) 

- 2s 2 (t + w)(16t 6 + 180t 5 M + 495* V + 672tV + 495* V + 180£m 5 + 16m 6 ) 

- s(t + m) 2 (7* 6 + 123t 5 M + 376tV + 465tV + 376tV + 123*m 5 + 7m 6 ) 

- 2tu(t + u) (t 2 + tu + u 2 ) (7t 4 + 48t 3 M + 76* V + 48*m 3 + 7m 4 )) 

+ (* + m)(s 6 (* + m)(2* 2 + Utu + 2m 2 ) + s 5 (4* 4 + 47t 3 M + 70* V + 47*m 3 + 4m 4 ) 

+ s 4 (t + m)(2* 4 + 6H 3 m + 75* V + 61*m 3 + 2m 4 ) 

+ s 3 *m(47* 4 + 132* 3 m + 190*V + 132*m 3 + 47m 4 ) 

+ s 2 *m(* + m)(25* 4 + 63t 3 M + 93* V + 63*m 3 + 25m 4 ) 

+ s*m(7* 6 + 38t 5 M + 78* 4 m 2 + 98t 3 M 3 + 78*V + 38*m 5 + 7m 6 ) 

+ 7*V(* + m)(* 2 + tu + m 2 ) 2 )] 

- 10a 2 Q s[36m 2 z ° - 8m^ 8 (16s + 17(* + u)) 

+ mf(212s 2 + 320s(* + u) + 153* 2 + 382*m + 153m 2 ) 

- m^ 4 (318s 3 + 346s 2 (* + u) - 2s(ll* 2 - lUtu + 11m 2 ) 

- (* + m)(67* 2 - 150*m + 67m 2 )) 

+ m^ 2 (440s 4 + 636s 3 (* + u) - s 2 (373t 2 + 362*m + 373m 2 ) 

- s(t + m)(915* 2 + 1034*m + 915m 2 ) 

- 337* 4 - 949* 3 m - 1180*V - 949*m 3 - 337m 4 ) 

- m^°(388s 5 + 982s 4 (* + u) + 2s 3 (128* 2 + 417*m + 128m 2 ) 

- s 2 (t + m)(1351* 2 + 1864*m + 1351m 2 ) 

- 2s(687* 4 + 2335t 3 M + 3216*V + 2335*m 3 + 687m 4 ) 

- 2(* + m)(187* 4 + 636* 3 m + 818* V + 636*m 3 + 187m 4 )) 
+ m|(188s 6 + 718s 5 (* + u) + 2s 4 (421* 2 + 946*m + 421m 2 ) 

- 2s 3 (t + m)(151* 2 + 45tM + 151m 2 ) 

- 2s 2 (703* 4 + 2496t 3 M + 3508*V + 2496*m 3 + 703m 4 ) 

- s(t + m)(1003* 4 + 3708* 3 m + 5034* V + 3708*m 3 + 1003m 4 ) 

- 217* 6 - 1323* 5 m - 3116tV - 4016* 3 m 3 - 3116* V - 1323tM 5 - 217m 6 ) 

- m|(46s 7 + 240s 6 (t + u) + 5s 5 (101t 2 + 220tM + 101m 2 ) 
+ s 4 (t + M)(339t 2 + 916tM + 339m 2 ) 
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- 2s 3 (196t 4 + 627t 3 M + 833tV + 627tu 3 + 196m 4 ) 

- s 2 (t + m)(77K 4 + 2854t 3 M + 3968tV + 2854tM 3 + 771m 4 ) 

- s(421t 6 + 2566t 5 M + 6173tV + 8040£ 3 m 3 + 6173tV + 2566£m 5 + 421m 6 ) 

- (t + u)(72t 6 + 524t 5 M + 1295tV + 1674* V + 1295tV + 524£m 5 + 72m 6 )) 
+ m z (As 8 + 30s 7 (t + u) + 2s 6 (53t 2 + lUtu + 53m 2 ) 

+ s 5 (t + m)(163£ 2 + 388£m + 163m 2 ) + s 4 (29t 4 + 237t 3 M + 424tV + 237tM 3 + 29m 4 ) 

- 2s 3 (t + M) 3 (108t 2 + 127tw + 108m 2 ) 

- s 2 (254t 6 + 1423t 5 M + 3385tV + 4416tV + 3385tV + 1423£m 5 + 254m 6 ) 

- s(t + M)(105t 6 + 682t 5 M + 1694t 4 M 2 + 2206tV + 1694t 2 M 4 + 682tM 5 + 105m 6 ) 

- (t + M) 2 (13t 6 + 126t 5 M + 331tV + 427tV + 331tV + 126£m 5 + 13m 6 )) 

- m 2 z (5s 7 (t + m) 2 + s 6 (t + M)(17t 2 + 38tM + 17m 2 ) 
+ s 5 (llt 4 + 67t 3 M + llOtV + 67tM 3 + 11m 4 ) 

- s 4 (t + u)(3t 2 + htu + 3m 2 ) (lit 2 + 2tu + 11m 2 ) 

- s 3 (69t 6 + 310t 5 M + 677tV + 866tV + 677tV + 310£m 5 + 69m 6 ) 

- s 2 {t + M)(51t 6 + 252t 5 M + 598tV + 776tV + 598tV + 252tM 5 + 51m 6 ) 

- s(t + M) 2 (15t 6 + 99t 5 M + 259tV + 336t 3 M 3 + 259tV + 99£m 5 + 15m 6 ) 

- (t + u) 3 (t 6 + 17t 5 M + 48t 4 M 2 + 62t 3 M 3 + 48t 2 M 4 + 17tM 5 + m 6 )) 

-(t + u)(s + t + u) 2 {s\t + u)(t 2 + tu + u 2 ) + s 3 (3t 4 + 6t 3 M + 8fV + 6tM 3 + 3m 4 ) 
+ s 2 (t + u) (3t 4 + 4t 3 M + 9tV + 4tM 3 + 3m 4 ) 

+ S {f + 4t 5 M + lltV + 14t 3 M 3 + lltV + 4tM 5 + M 6 ) 

+ tu{t + uf{t 2 + tu + u 2 ))}}, (A.14) 



da 
~dt 



(qq 



QQ psf] Z) 



9Ms 3 (2m| -t -u) 



-[2m| - 2m 2 z (s + t + u) + t 2 + m 2 ], 



(A.15) 



da / _ 



QQ 6 S\ 



3q(8) 



2^2 



-2najg 



27M 5 m|s 3 t% 2 (2m| 



u 



x {{v 2 + a 2 )m 2 z s(2m 



[m z {t 2 



u 



2m\ 



t + u){t 2 + m 2 ) 



+ m 2 z {s 2 {t + m) 2 + 2s(t + u){t 2 + tu + m 2 ) + (t 2 + M 2 )(t 2 + 3tM + m 2 )) 

- tM(2s 3 + 4s 2 (t + u) + s(3t 2 + Atu + 3m 2 ) + (t + u){t 2 + m 2 ))] 

— Aa g aQM 2 m 2 z stu(2m 2 z — t — u)(t + u)[2m z — 3m|(s + t 



u 



(s + t + u) 



+ 2a 2 Q MH 2 u 2 [2m z (3s - t - u) - m 2 z (Qs 2 + As(t + u) - 2(t + u) 
- (s + t + M)(t 2 + M 2 )]}, 



(A.16) 
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da 
~dl 

x {8m|(s 



(gg -> gg 



(8) 



— 167ra: 2 g 2 a 
27M 3 m|s 3 (2m| 



Q 



i-u) - 16m|[s 2 - (£ + m) ; 



+ 2m|[4s 3 + 4s 2 (t + «) - s(3t 2 + lOfu + 3m 2 ) - (t + u)(7t 2 + Wtu + 7m 2 )] 
- 2m 2 z [s 2 (t - u) 2 - 2s(t + m) 3 - (t + u) 2 (3t 2 + 2tu + 3m 2 )] 

2/ 4-1 



(s + t + u)(t + u) 2 (t 2 + m 2 )}, 



(A.17) 



da / 



gg 



3p 



(8) 
J 



leTra 2 ^ 



9M 3 s 3 (2m| 
x {16m! - 8m|[8s + 5(f + u)} + 2m|[16s 2 

+ 2m|[8s 3 - 4s 2 (i + u) - 23s(t + m) 2 - (t + m)(13^ + 14fu + 13m 
+ (t + u)[8s 2 {t + m) + 4s(3t 2 + Atu + 3m 2 ) + 7(t + m)(£ 2 + m 2 )]}, 



48s(t + u) + 23t 2 

(-2 



38tM + 23m' 

2\1 



(A.18) 



da 
~dt 



(du - gg [ 3 s{ 8) ] w-) = § (ud -> gg [ 3 s{ 8) ] w+) 



27M 3 s 2 t 2 M 2 



X {m^(t 2 + M 2 ) - m^(s + i + M)(t 2 + m 2 ) + tM[2s 2 + 2s(t + u) + t 2 + m 2 ]}. (A.19) 
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ee -> ee C D + X at TESLA in ee mode, C = J/\|/, %, 
p T [GeV] p T [GeV] 



5 10 15 20 5 10 15 20 




5 10 15 20 5 10 15 20 

P T [GeV] p T [GeV] 

(a) 

Figure 1: (a) pr distributions da/dpx (in fb/GeV) and (b) yc distributions da/dyc (in fb) 
of e + e~ -> eVCT + I, where C = J/ij), Xc j and D = 7, Z, //, at TESLA in the e+e" 
mode. It is summed over C = Xco,Xci,Xc2 and D = W + ,W~. In each figure, the CSM 
(dotted lines) and NRQCD (solid lines) predictions for C = J/ip and the CSM (short- 
dashed lines) and NRQCD (medium-dashed lines) ones for C = Xcj are shown separately. 
Superposed short-dashed and medium-dashed lines appear as long-dashed lines. 
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ee — > ee C D + X at TESLA in ee mode, C = JA|/, % 




e e — > e e C D + X at TESLA in yy mode, C = J/\|/, % cj 
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Figure 2: Same as in Figs. |l|(a) and (b), but for e + e — > e + e CD + X at TESLA in the 
77 mode. 
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e e -> e e C D + X at TESLA in yy mode, C = J/\\f, 




pe *€-> e ± €C D + X at THERA, C = J/\|/, %, 



P T [GeV] p T [GeV] 
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P T [GeV] p T [GeV] 



(a) 

Figure 3: Same as in Figs. [l|(a) and (b), but for pe ± e ± CD + X at THERA in the 
laboratory frame. 
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pe e *€C D + X at THERA, C = JA|/, X, 
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pp — > C D + X at the Tevatron, C = JA|/, X, 





40 



pp -> CD + X at the LHC, C = JA|/, %, 




